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Abstract 

At the short times, the enstrophy 17 of a two-dimensional flow, generated by a random 
Gaussian initial condition decays as ^}{t) oc t''^ with 7 ~ 0.7. After that, the flow 
undergoes transition to a different state characterized by the magnitude of the decay 
exponent 7 0.4 (Yakhot, Wanderer, Phys.Rev.Lett.93, 154502 (2004)). It is shown 
that the very existence of this transition and various characteristics of evolving flow 
crucially depend upon phase correlation between the large- scale modes containing only 
a few percent of total enstrophy. 
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Two-dimensional (2D) turbulence is a simplified system mimicking such important phe- 
nomena, as ubiquitous vortex generation, in atmospheric flows. The vortex formation as a 
dominant process of 2D hydrodynamics has been recognized by Onsager more than half a 
century ago^ and, thanks to obsession of TV stations with the weather forcast, the relevance 
of spontanious emergence of strong vortices out of quiescent background for the world we live 
in, can be verified at any time of the day. Many features of the vortex formation in decaying 
2D turbulence have been undestood due to numerical and theoretical works of Mc Williams 
and collaborators^"^, Benzi et. al.^ and many others. Still, many new features have recently 
been discovered and some important questions remained unanswered. 

The problem of decay of two-dimensional (2D) turbulence considered in our recent papcr^ 
and the one, wc are interested in here, is formulated in a following way. Consider a time 
evolution of an initial velocity field, u{x, y, 0) = Uq defined on a two-dimensional square such 
that —L < x,y < L. The field uq is a Gaussian random noise supported in the Fourier 
space in the vicinity of A; = ko{t = 0) = O(^) with = 0(1) and N = const = 0(1). At 
the initial instant, t = 0, the enstrophy is Q = o;^ = cUq = (V x u)^ = 0(1). The square 
is large meaning that we are interested in the limit u ^ and ko{t = 0)L oo. Still, 
the box is finite so we will be able to study both short time, when ko{t)L ^ 1, and long 
time asymptotics when ko{t)L 1. In all our simulations the initial kinetic energy of the 
flow K{t — 0) — J v'^{x.,t — 0)dx = i and the fluid density p — 1. This set up with 
the basically structureless initial fleld is ideal for investigation of the details of the structure 
emergence in the course of the flow evolution. 

According to a recent theory®, in an inflnite system (L oo), the universal asymptotic 
law of enstrophy decay is Q oc t~'^ with 7 = | which is close to the short -time results of 
simulation by Chasnov^ (7 f« 0.7 — 0.8) and our own''' 7 Ri 0.7. At the later times, after 
the large-scale (small-wave-vector) modes are somewhat populated, the system undergoes 
transition to a gas or liquid of well separated vortices and, simultaniously, the magnitude of 
the exponent 7 crosses over^ to 7 ~ 0.4, well-known from studies of evolution of an initially 
created ensemble of point vortices^"^. The theory developed in Ref.[7] gave for the number 
of vortices in the system N{t) oc with ^ = 4/5 and 7 = 2/5 very close to numerical 
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data. This theory was based on the work by Carnevale et aF which led to the relation 
7 = ^/2 with an undetermined magnitude of exponent ^. Although a tentative correlation 
between the cross-over to 7 0.4 with population of the large-scale modes was mentioned in 
Ref . [7] , the details of the process remained obscure. All theories treated the vortex merger 
process locally as a binary "chemical reaction" 2n ^ n with the constant peak vorticity 
before and after collision. The role of the large-scale patterns on scales / = 0{L) in this 
process was completely ignored. Below we report a surprising discovery: the very existence 
of transition to the long-time asymptotic (7 = 0.4) of the enstrophy decay and the energy 
of evolving vortices crucially depend upon phase correlation (coherence) of the large-scale 
modes containing only a tiny fraction of the total enstrophy O. 

The Navier-Stokes equations with the 0(z/4V^v) hyper- viscous dissipation terms were 
simulated using a pseudo-spectral method. The initial random field was Gaussian with 
... .e.™. - .4 ... . . ^ a. . . 3 . . 

Ref. [9] . The parameters of the simulations are given in a table of Ref. [7] . 

On Fig. 1 the time evolution of normalized enstrophy Q{t)/Q{t = 0) is presented. The 
curve A is the one obtained from the 512^ simulations reported in Ref. [7]. The cross-over from 
from the short-time magnitude of the decay exponent 7 ~ 0.7 to the long -time asymptotics 
7 0.4 happening at the dimensionless time T — tcoo 120 is clearly seen (curve A). To 
investigate the role of the large-scale modes in this phenomenon we, at the transition time 
T ^ 120, randomized the phases of the first few modes {k < 5 ) (curve R) and let the flow 
evolve. It is important to stress that this procedure influenced neither total enstrophy nor 
energy in the system. Still, the observed response of the flow to this perturbation was quite 
dramatic: the randomization of phases of the modes, containing only a few percent of the 
total enstrophy, prevented transition from happening. Then, following the discussion with 
J.McWilliams, instead of randomizing, we at the same instant T = 120 inverted the signs 
of the same modes in the unperturbed run A, i.e. made a transformation Vi{k) —Vi{k) 
with k < 5. The effect of this procedure on the enstrophy evolution can be seen on Fig. la 
(curve N): the curves R and N are almost indistinguishable. After some time {T 470 ) the 
perturbed flow (R) developed the large-scale coherent motions and a tendency to crossover 
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to the long-time asymptotics (7 ~ 0.4) could be detected. At this instant , again, the phases 
of the modes with A; < 5 have been randomized. The result of this perturbation is shown on 
the curve RR of Fig. la. We can see that the second randomization too, forced the enstrophy 
to decay much faster than D, oc observed in the dynamically unperturbed run A. At this 
point we can make a definitive statement: the phase correlation , dynamically established 
as a result of the large-scale evolution, is crucial for the very existence of transition to the 
long-time enstrophy decay regime. 

If the previous experiment established strong coupling between large {k < 5) and small ( 
enstrophy containing) scales , then the next natural question to ask is: what is the minimal 
length-scale of the energy-containig structures influencing such small-scale phenomenon as 
the enstrophy decay process? To answer this question, we broadened our search. The phase 
randomization of the modes v{k) with k < 2 did not show any effect on the decay process. 
The result of randomizing the phases of the modes with k < A containing at T 120 less 
than five percent of total enstrophy Jl, is shown on Fig. lb: the transition to the long -time 
behavior is delayed. 

The time-dependence of the fourth and sixth-order moments are shown on Fig. 2. As 
we see, the initially gaussian random flow {t — 0) develops strongly non-gaussian tails char- 
acteristic of the small-scale vortex formation also observed in Ref.[10]. It is interesting that 
at the long times, the normalized high-order moments of the unperturbed vorticity field 
(A) approach a close to steady state. Thus, in in accord with Ref.[ll], we can write an 
exact equation for probability density P(X) = P(— ^) in terms of conditional expecta- 
tion value of vorticity dissipation rate Q{X) for a given magnitude of X. A simple model 
Q{X) oc (1 -I- kX'^) proposed and numerically verified in Refs. [11] gives: 

PiX) = r (1) 

where C is a normalization constant. After the transition (T > 120), the maximum value 
of vorticity in the evolving unperturbed field cumax/^^o ~ 60 ~ const and the expression (1) 
is valid in the interval Xmax < ^^rnax/^^rms{t). For X > X„iax, the probability density is cut 
off, i.e. we can set P{X) fti 0. It is due to the cut-off X^nax, the probability density (1) is 
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consistent with the finite moments of vorticity field. Since oOrms = the magnitude 

of X^ax oc does not appreciably varies during the time of our simulation and can be set 
Xmax ~ const. The data (curves A) reported on Fig. 2 are reasonably well fitted by k 1 
and Xjnax ~ 30. As k ^ 0, the relation (1) approaches a Gaussian and, since the evolution 
of the normalized moments is slow, the curves A of Fig. 2 can be represented by (1) with the 
time-dependent parameter K{t). As we see from Fig. 2, during the time of simulation, the 
sixth-order moments of perturbed runs R and do not reach steady state meaning that the 
strong vortices decay with the exponent 7^ < 0.7. 

To make the discussion easier we will be dealing with three fields: 1. unperturbed (A); 2. 
randomized (R); and 3. sign-reversed (N). The time-evolution of vorticity field in all three 
runs (A; R N) is shown on Figs. 3-5. On Figs.3(A;R;N;) the vorticity fields at the transition 
time T 120 right after randomization (Fig.3R) and sign reversal (Fig. 3N) are compared 
with the original unperturbed field (Fig. 3A). In all three plates we can identify the small- 
scale strong vortices which have not been affected by the randomization or sign-reversal. 
The patterns between these structures are quite different, though. The time-evolution of 
these fields presented on Figs.4-5 for the times T k, 480 and T k, 975, respectively, leads to 
substantially different patterns: the larger number of strong vortices (also see below) can be 
seen in the field A, while the evolution of the fields R and N results in the smaller number 
of vorticies. While during this evolution, the maximum value of vorticity stayed more or 
less unchanged {oJmax ~ 65 — 60), the number of strongest vortices having the peak value of 
vorticity uJo ~ O.ScJmaa; (not shown here due to lack of space) in the unperturbed field A was 
much larger (factor 1.5-4, depending on time T) than that in the corresponding perturbed 
runs R and A^. 

In Fig. 6 we show the time-evolution of maximum velocity v^nax in three cases as 
a function of time. We can sec that after some transient, the velocity in the unperturbed 
field A reaches the values v„iax ~ 3.5 — 4., while in the time interval 120 < T < 500 the 
velocity in perturbed fields R and N is smaller: v^ax ~ 2.7 — 3.. It is interesting that by the 
time T fa 600, the maximum velocity of the field R is more or less recovered while that in 
a stronger perturbed field N remained substantially smaller than that in the A-field during 
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the remaining time of the simulation. 

The strong correlation between small and large-scale structures in the small-scale forced 2D 
turbulence was discovered in Refs. [12]-[14]. Two different flow regimes were identified: at 
the short times, when the time-dependent integral scale Ij <^ L, the inverse cascade led 
to generation of two -dimensional turbulence characterized by the Kolmogorov spectrum 
and close to gaussian magnitudes of the even-order moments velocity differences. The most 
interesting effect happened at the later times when the modes with A; ~ 1/L became pop- 
ulated: very strong vortices, were simultaneously formed at the forcing scale If ^ 1/kf. It 
was demonstrated^^ that in the case of forced turbulence, the vortices were created at the 
centers of the large-scale, slowly -varying patterns where velocity v fa 0. Thus, the Bose 
condensate served as an almost steady well-organized matrix facilitating the vortex merger 
process. The possible generality of this phenomenon was demonstrated by the experiments 
of Paret and Tabeling^^ in their investigation of time-evolution of the two-dimensional flow 
of mercury generated by magnetic field. 

It is important to stress that, unlike in the forced flow considered in Refs. [12]-[14] where 
vorticity uj^ax of the strongest vortices continuously grew with time, in decaying turbulence 
studied here, oJmaxif) ~ const ~ 60. If this feature persists for a very long time, eventually 
only two strong vorticis of the radius a^o ~ [K / {poj'^g^^)]^ will be left in agreement with 
Ref.[l]. For the flow considered in this paper aoo ~ 0.1. Thus, the final vortex occupies a 
tiny fraction (k, 10~^ — 10~^) of the flow. The flow velocity in this configuration will reach 
V ~ doo^max ~ \J'^max [K / p] 4 ^ 5 — 10 which is much larger than that in the initial random 
state. Thus, the 2D turbulence decay leads to formation of localized high-energy structures 
out of the initially quiescent background. 

To conclude: the main result of the present study is: we have showed that the cross-over 
between the two asymptotics of enstrophy decay in 2D turbulence strongly depends upon 
fine details of the large-scale flow features containing very small fraction of total enstrophy 
n. The dynamic cause of this transition is not yet clear. However, we can conclude that fine, 
dynamically evolved, phase-correlation of the large-scale modes leads to a slower enstrophy 
decay and formation of more energetic flow structures. This points to a strong coupling 
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between the small-scale and large-scale (condensates) structures in 2D turbulence. At the 
present time, we do not have a theory leading to this coupling and do not fully appreciate 
it's consequences. Based on previous examples of coherent-random states interaction (super- 
fluidity, superconductivity and many others) , the theory of the effect observed in this paper 
may be a substantial theoretical challenge. The possibility of this coupling in forced 2D 
turbulence was discussed by Polyakov^^ in his conformal theory of two-dimensional turbu- 
lence. Since the geometry-dependent large-scale patterns cannot be universal, the question 
of universality of the transition is extremely interesting. The possible role of the inter-scale 
interaction, observed in this work, may be of great interest for meteorological applications. 
We are grateful to J. McWiUiams and A.Polyakov for valuable comments and suggestions. 
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(a) Time evolution of total enstrophy Q. A. Unperturbed field. R. ef- 
fect of randomization of the modes with fc < 5; N. n{t) at T > 120 
with v{k) — > —v{k) {k < 5); RR. evolution after second randomization 
at T w 470. 
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(b) A. The same as on la; R4: Evolution after randomization of the 
modes with k < 4. 



Figure 1: Enstrophy vs Time 
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Figure 2: Time-evolution of vorticity moments for the three runs 
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(a) A; Lu{x,y); T=120. 



(b) R;uj{x,y); T=120 



(c) N; c<j(a;,y)|/c| < 5 



Figure 3: Vorticity fields A, R, and N at the transition time 
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Figure 6: Time-evolution of maximum velocity for three runs: Curve ^4: Vmaxif)] R'- Vmax + ^ 
vs. t; N: v^ax + 2 vs t; 
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